Critical adsorption and Casimir torque in wedges and at ridges 
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Geometrical structures of confining surfaces profoundly influence the adsorption of fluids upon 
approaching a critical point T c in their bulk phase diagram, i.e., for t = (T—T c ) /T c — > ±0. Guided by 
general scaling considerations, we calculate, within mean-field theory, the temperature dependence 
of the order parameter profile in a wedge with opening angle 7 < 7r and close to a ridge (7 > n) 
for T ^ T c and in the presence of surface fields. For a suitably defined reduced excess adsorption 
r±(7,t — > ±0) ~ T±( r y)\t\ l3 ~ 21 ' we compute the universal amplitudes T±( r y), which diverge as 
r±(7 — * 0) ~ I/7 for small opening angles, vary linearly close to 7 = n for 7 < 7r, and increase 
exponentially for 7 — > 2tt. There is evidence that, within mean-field theory, the ratio r+(7)/r_(7) 
is independent of 7. We also discuss the critical Casimir torque acting on the sides of the wedge as 
a function of the opening angle and temperature. 

PACS numbers: 64.60.Fr, 68.35.Rh, 68.43.Fg, 61.20.-p 
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I. INTRODUCTION 

Boundaries induce deviations of the local structural 
properties of condensed matter from their correspond- 
ing bulk values. Typically the width of such boundary 
layers is proportional to the bulk correlation length £. 
Near a continuous bulk phase transition at a tempera- 
ture T = T c , the correlation length diverges according 
to a power law £ ± (t = (T — T c )/T c -> ±0) = ^\t\~ u 
with a universal bulk exponent v and nonuniversal am- 
plitudes £q whose ratio £o~/£cT 1S universal, too. For 
the three-dimensional Ising universality class v ~ 0.63 
and £o~/£o~ -2 0, whereas within mean-field theory, 
which is valid for spatial dimensions d > 4, v = 1/2 and 
£o~/£o~ = V^2- The nonuniversal amplitudes are typi- 
cally in the order of the range of the interaction potential 
of the ordering degrees of freedom, i.e., a few A. The cor- 
relation length is defined as the scale of the exponentially 
decaying two-point correlation function ~ e~ r ^. 

The local critical properties near planar confining sur- 
faces have been studied theoretically and experimentally 
in detail. It has turned out, that each bulk universal- 
ity class splits up into three possible surface universality 
classes denoted as ordinary, special, and normal transi- 
tions. Each of them gives rise to distinct surface critical 
phenomena. The type of boundary conditions determines 
which surface universality class a given system belongs 
to. The ordinary transition requires the absence of sur- 
face fields, the special transition is characterized by the 
absence of a surface field, too, but also by suitably tuned, 
enhanced surface couplings between the ordering degrees 
of freedom, and the normal transition represents systems 
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whose surfaces are exposed to surface fields. As indicated 
by the name, the latter one is the generic case. For exam- 
ple the normal transition applies to one-component fluids 
near their liquid- vapor critical point, or to binary liquid 
mixtures near their demixing critical point. In both cases 
the substrate potential of the confining container walls 
and the missing interactions due to the fact that the fluid 
particles cannot penetrate the substrate give rise to the 
effective surface fields acting on the corresponding or- 
der parameter (density or concentration difference). As 
a consequence, in this case even above T c there is a non- 
vanishing order parameter profile, which decays to zero 
for T >T C upon approaching the bulk, i.e., for increasing 
normal distance z — > 00 from the surface at z — 0. For 
T < T c this profile attains the nonzero value of the bulk 
order parameter for z — ► 00. Following early calculations 
on magnetic systems by K. Binder and P. C. Hohenberg 
, this so-called critical adsorption was first examined in 
detail by Fisher and de Gennes and has since been an- 
alyzed for many systems both theoretically 0, 0, 01 and 
experimentally [3, and a fair agreement between theory 
and experiment has been found 0, |j| . 

Surface order above the bulk critical temperature can 
also be due to spontaneous symmetry breaking caused 
by enhanced surface couplings between the ordering de- 
grees of freedom, e.g. spins in magnetic systems. The 
transition to this state from the disordered state is usu- 
ally denoted as the extraordinary transition. Bray and 
Moore [Tol ] predicted an equivalence between the normal 
and the extraordinary transitions that was later proved 
by Burkhardt and Diehl In contrast to the or- 

dinary transition in magnetic systems, the correspond- 
ing extraordinary transition has been investigated to a 
lesser extent. Extending and impro ving earlier work 
011113,0, Diehl and Smock [j]f have carried 
out a field-theoretic renormalization-group calculation in 
4 — e dimensions for the extraordinary transition in semi- 
infinite systems belonging to the Ising universality class 
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[lTj computing the order parameter profile to one-loop 
order. The results are in fair agreement with those ob- 
tained by Monte Carlo simulations Q. 

These studies are devoted to the case of planar sur- 
faces. Only very carefully treated solid surfaces are atom- 
ically flat; generically, however, they exhibit corruga- 
tions. Besides these random deviations from the flat to- 
pography, there is nowadays an abundance of experimen- 
tal techniques 0| that allow one to endow solid surfaces 
with precise lateral geometrical structures, ranging from 
the \xm scale down to tens of nm. Inter alia, such solid 
surfaces are used within the context of microfluidic de- 
vices fl9| in order to guide fluids along such structures. 
The fluids perfectly fill this laterally structured environ- 
ment and thus fully exhibit the strong structural changes 
associated with that. In this context it is of interest what 
kind of structures appear in the fluid if it is brought, say, 
close to a demixing transition and is exposed to geomet- 
rically structured substrates. If the characteristic sizes 
of the lateral structures are comparable with the correla- 
tion length, which in fluids reaches typically up to a few 
thousand A close to T c , one can expect a strong influence 
on the aforementioned critical adsorption phenomena. 

Theoretically this raises the issue of how the local crit- 
ical properties depend on the shape of the boundaries. 
A typical manmade structure is a series of grooves with 
various shapes of the cross section, e.g. wcdgclike. The 
first step in their investigation is the study of a single 
wedge, which already shows new features and gives new 
insight into the influence of geometry on critical behav- 
ior. The ordinary transition of the isotropic N- vector 
model at an edge has been investigated by Cardy |2(j 
within the framework of mean-field theory, the renor- 
malization group, and e-expansion. Subsequently the 
two-dimensional wedge (corner geometry) was studied by 
exact calculations, mainly for Ising models |2l|. and by 
conformal mapping at the bulk critical temperature [23] • 
New edge and corner exponents were found that depend 
on the opening angle 7 of the wedge. More recently sim- 
ilar findings of Monte Carlo simulations of the ordinary 
transition of the three-dimensional Ising model [23T ] have 
been reported, in accordance with high temperature se- 
ries expansions |24j ■ The angle dependence of the critical 
wedge exponents can be attributed to the fact that the 
wedge geometry lacks a length scale and thus is invari- 
ant under rescaling. The opening angle is therefore a 
marginal variable in a renormalization transformation, 
and may enter into the expressions for the exponents. 
The same will happen for all other scale-invariant shapes 
of the boundaries |25j . However, for a given opening 
angle, the values of the critical exponents arc expected 
to be universal and independent of microscopic details. 
According to recent Monte Carlo simulations of three- 
dimensional Ising models with edges and corners [2(|, 
angle dependent critical exponents are observed at the 
ordinary transition, whereas the surface transition seems 
to be nonunivcrsal. The critical exponents in this latter 
case appear to depend on the strengths of the local cou- 



plings, in analogy with exact results obtained for the two- 
dimensional Ising model with defect lines [27|]. Critical 
adsorption has also been studied in general dimensions d 
on curved surfaces psj . 

Since critical adsorption changes the composition of 
a binary liquid mixture close to its surface, mechanical 
properties, such as the local viscosity and the mutual 
diffusion coefficient will also vary in space. So we expect 
various phenomena, such as flow properties in porous me- 
dia, the spreading properties of droplets, surface chemical 
reactions, the permeability of membranes etc. p9j . to be 
influenced significantly by critical adsorption. 

Critical adsorption in a wedge has been studied by 
Hankc ct. al. ,30]. Within mean- field theory the order 
parameter profile was calculated exactly at the critical 
point. Through an interpolation scheme between exact 
results in d = 2 and the mean-field results correspond- 
ing to d = 4, angle dependent critical exponents of the 
order parameter and those governing the decay of two- 
point correlation functions were obtained for d = 3 as 
well. The present study extends these investigations into 
various directions. 

First, we analyze the temperature dependence off crit- 
ically for critical adsorption in a wedge. General scaling 
properties for the order parameter profile and the excess 
adsorption are discussed and the corresponding scaling 
functions are calculated within mean-field theory. This 
analysis is carried out above and below T c , including a 
systematic study of the dependence on the opening angle 
7 of the wedge. This covers also the case 7 > ir describing 
critical adsorption near a ridge. 

Second, as a new feature, we study the Casimir torque 
acting on the sides of the wedge or ridge. The confine- 
ment modifies the fluctuation spectrum of the critical 
fluctuations and the order parameter profiles. This leads 
to a dependence of the free energy of the critical medium 
on the shape of and the distance between the confining 
boundaries, which results in an effective force acting on 
them. Thus the physical origin of this force, originally 
predicted by Fisher and de Gennes for two parallel 
plates immersed into a binary liquid mixture near its con- 
tinuous demixing transition, is analogous to the Casimir 
force acting on conducting plates in vacuum due to the 
confinement of quantum mechanical vacuum fluctuations 
of the electromagnetic field [31). The Casimir force is 
governed by universal scaling functions |32T | and is super- 
imposed on the noncritical background forces, which in 
the case of fluids are given by dispersion forces. Recent 
experiments |33l ] have confirmed corres pon ding theoret- 
ical predictions for the plate geometry j34j. For curved 
surfaces the critical Casimir force plays an important role 
in the flocculation of colloidal particles suspended in a 
solvent undergoing a continuous phase transition |35j . In 
the present context the free energy of the critical medium 
depends on the opening angle 7; its derivative with re- 
spect to 7 amounts to the critical Casimir torque acting 
on the sides of the wedge or ridge. If the substrate form- 
ing the wedge or ridge is composed of soft materials like, 
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e.g., membranes, this critical Casimir torque is expected 
to give rise to elastic deformations. It might also be ex- 
perimentally accessible by suitable force microscopy with 
moveable sidewalls of a wedgelike structure. 

The paper is organized as follows. In Sec. II we dis- 
cuss the general scaling properties of the order parameter 
profiles and the excess adsorption. The scaling functions 
are determined within mean-field theory and analyzed in 
detail in Sec. III. Section IV focuses on the free energy of 
the confined fluid and the critical Casimir torque result- 
ing from its angle dependence. Section V summarizes our 
findings. In the Appendix we discuss how the excess ad- 
sorption in a wedge or at a ridge decomposes into surface 
and line contributions with two possible experimental re- 
alizations. 



II. GENERAL SCALING PROPERTIES OF 
ORDER PARAMETER PROFILES AND EXCESS 
ADSORPTION 

Since fluids can fill a container of arbitrary shape, in 
the present context of critical systems exposed to sub- 
strates shaped as wedges we consider fluids close to their 
bulk critical point T c . This can be either a liquid- vapor 
critical point or a demixing critical point in the case of 
binary liquid mixtures. The interaction of the container 
walls with the fluid particles results in a spatial variation 
of the number densities close to the boundaries. The de- 
viation of the density of the fluid, or of the concentration 
of one of its two components in the case of binary liq- 
uid mixtures, from the corresponding bulk value at T c is 
chosen as the local order parameter describing the phase 
transition. 

The order parameter profile m±^ 2 (£,t) near a planar 
interface and close to the critical temperature T c takes 
the following scaling form H H Q E : 

m^ /2 (U)=a\tfP^ /2 ((/Z ± ), t={T-T c )/T c , (1) 

for distances £ > a perpendicular to the interface and suf- 
ficiently large compared to a typical microscopic length a. 
£±(t — ► 0) = ^o = |i| -1 ' is the bulk correlation length above 
(+) or below (— ) T c , (3 and v are the standard bulk crit- 
ical exponents, and t is the reduced temperature. The 
scaling functions P±^ 2 (C± = C/£±) are universal once 
the nonuniversal bulk amplitudes a and £q are fixed, 
where a is the amplitude of the bulk order parameter 
m!° /2 (C = oo, t ->■ 0") = a\tf = m b (t). With the pref- 
actors £^ fixed as those of the true correlation length de- 
fined by the exponential decay of the bulk two-point cor- 
relation function in real space, one finds P^°^ 2 (oo) = I , 

/2 (oo) = 0, P°°/ 2 (Q_ -» oo) - 1 ~ e-f-, P+ /2 (C+ 
oo) ~ e~<+ and P^ /2 (C± -► 0) = c±(g Pl " 0, so that 

m ac/2 (C,t = 0) = ac ± (C/^r^- (2) 



Any other choice for the definition of the correlation 
length leads to a redefinition of the scaling functions 
poo/2 guc j 1 j^at all observable quantities remain un- 
changed. This underscores that the scaling functions are 
universal, but that their form depends on the definition 
of the correlation length. The amplitudes of the scaling 
functions arc fixed by the requirement P^°^ 2 (oo) = f. 
Accordingly the numbers c± arc universal surface ampli- 
tudes which are definition-dependent @. 

Close to T c the total enrichment at the interface of, 
say, the A particles as compared to the B particles of a 
binary liquid mixture is given by the excess adsorption, 
which is an experimentally accessible integral quantity. 
For a planar surface, one has 

/>oo 

T™ /2 (t)= {m™ /2 ((,t)-m™ /2 (( = ^t)}dt (3) 
Jo 

The scaling behavior of T±^ 2 (t) has been discussed in 
Ref. Q in detail. The £ integration can be split into the 
intervals ( > a and < £ < <r, and for large £ the order 
parameter profile m± ^ 2 (C, t) can be replaced by Eq. QJ, 
which gives 

Yf 2 {t) = / CT [m^ /2 (C,t) - m^ /2 (C = oo,t)]dC 
Jo 

+ a£\tf-" / [P^ /2 (C±)-^ /2 (C ± =oo)]dC±. 

(4) 

Upon approaching T c the first integral remains finite and 
yields a nonuniversal constant, which is subdominant to 
the diverging second term. The second integral leads to 
the well known power-law sing ularity ofr± /2 (i^0) for 
d < 4: 

T^ 2 (t^0)=a^9± 1 -^, rf<4, (5) 

where the numbers g± are universal with their values 
depending on the definitions of the bulk order parameter 
and the correlation length. In d = 4 one finds upon 
inserting the mean-field expression for p±l 2 (see Ref. 

[l6j ') into Eq. Q that T±^ 2 (t) diverges logarithmically. 
This result can be reconciled with Eq. (JSJ by noting that 
Y°?l 2 actually needs additive renormalization leading to 

a 

Tt ,2 {t^Q) = aitg J^~ V - 1 , d < 4. (6) 



In the wedge or ridge geometry (Fig. the order pa- 
rameter profile depends on the radial distance r from the 
apex, on the polar angle 9, and on the opening angle 7, so 
that the variation of the profile is two-dimensional with 
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corresponding generalized scaling functions: 

m±(r,i;7)=a|t|0p±(r/fcfc,0;7). (7) 
As before, the scaling functions are universal once the 




ridge 

FIG. 1: A wedge and a ridge with opening angle 7 exposed 
to a binary liquid mixture. The system is translationally in- 
variant in the subspace parallel to the edge. Within a plane 
orthogonal to the edge the polar coordinates are r and 9. The 
linear extensions of the confining surfaces are sm and s±. 

nonuniversal bulk amplitudes a and £^ are fixed, where 
a is the amplitude of the bulk order parameter and 
the amplitudes £ are prefactors of the true correlation 
length as in the case of the infinite planar wall. One 
finds [H3 with r± = r/£± that P_(Y_ = 00, 0) = 1, 
P+(r+ = oo,0) = 0, P_(r_ — > 00, 0) — 1 ~ e -sin(0)r_ 
(0 < 7/2), P+(r+ 00,8) ~ e - sin W r + (0 < 7/2), 
P±(r± 0,0) = c ± (0,7)r ± /3/ ' y , and 

P±(r ± ,0^O) = c ± C^ /1 ' 

= c±(r± sin0)-^, 

where C± = r±sin0 (see Fig. The amplitudes of 

the scaling functions are again fixed by the requirement 
P_(r = oo,0) = 1. The numbers c± are the universal 
surface amplitudes of the scaling function of the infinite 
planar wall [Eq. 0], and c±(9, 7) are universal func- 
tions. Both c± and c± depend on the definition of the 
correlation length. The scaling functions also reflect the 



symmetry of present the geometry: 

P±(r ± ,0; 7 ) = P±(r ± ,7-0; 7 ). (9) 

We define the excess adsorption for this geometry con- 
fined by surfaces of linear extensions sj_ in the plane per- 
pendicular to the edge and sn in the translationally in- 
variant directions (see Fig. ^) as 

f±(a±,*||,*;7)= / d d r[m±{r,t ;l )-m b {t)}, (10) 
Jv 

where the integral is taken over a macroscopic volume V 
occupied by the liquid. According to the Appendix this 
excess adsorption decomposes into a surface contribution 
that scales with the actual surface area of the confining 
walls (sj_S|| d 2 ^) and a line contribution that scales with 

the extension in the invariant directions (sj| d , i.e., a 
line in d = 3) j^: 

f ± ( S± , S ||,i; 7 ) = Tf(t) S± s\^ + Tf(t, 7 )s\^ 

+0(sl 1 ). (II) 

While the surface term is determined solely by the or- 
der parameter profile of a semi-infinite binary mixture 
exposed to a flat substrate, the line contribution is the 
specific contribution arising from the change of the order 
parameter profile caused by the edge. Due to the sym- 
metry of the system one can determine the amplitudes 
and for the wedge and the ridge explicitly by con- 
sidering only one half of the wedge or the ridge [see Eq. 
©], and by suitably subtracting and adding the order 
parameter profile m^ 2 (£ = rsin0,i) of a fluid in con- 
tact with an infinite planar wall in the integrand of Eq. 
(|10|l (compare Eqs. 10)-© and the Appendix): 

rf(t) = r± /2 (t) (12) 

and 

r±(t )7 ) = 2 / d9 drr[m±(r,0,t;i) 
Jo Jo 

-m^ /2 (C(r,0),i)]- (13) 
Based on Eq. Q T^(t, 7) takes on the scaling form 

r?(t, 1 ) = a£ 2 \tf- 2 »T ± ( 7 ) (14) 

with the universal amplitude functions 

/•7/2 poo 

r±( 7 ) - 2 / d0 dr ± r ± [P ± (r ± ,9;j) 
Jo Jo 

-P± /2 ((±(r±,e))} (15) 

where (±(r±, 0) = r± sin0. We note that the integral in 
Eq. (|15f) is finite for d = 4, i.e., T± does not carry a factor 
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proportional to as does (compare Eq. JSJl). 

As one can see from Eqs. © and 114fl . the subdominant 
line contribution to the excess adsorption carries a more 
singular temperature dependence than the leading planar 
surface term. Thus the scaling properties of the order 
parameter profile completely fix the functional form of 
the excess critical adsorption up to the dependence of the 
universal amplitudes T±(7) on the opening angle. Since 
this dependence cannot be inferred from general scaling 
arguments, it must be determined explicitly. This will 
be carried out within mean-field theory in the following 
section. This is possible because as stated above r±(7) 
is finite for d = 4. 



III. SCALING FUNCTIONS WITHIN 
MEAN-FIELD THEORY 

A. Order parameter profiles 

The standard Ginzburg-Landau Hamiltonian for de- 
scribing critical phenomena in confined geometries is 

Si 



(16) 



with a scalar order parameter field </>(»■•), supplemented 
by the boundary condition <fi = +°° a t the surfaces of 
the wedge (ridge) corresponding to the critical adsorp- 
tion fixed point . The parameter r is proportional to 
the reduced temperature t, u is the coupling constant, 
and the integration runs over the volume V w m of the 
wedge (ridge) (see Fig. Q). Within mean-field theory 
r = t/(&) 2 for T > T c , and r = -§|f|/(£-) 2 for T < T c 
with = y/2. 

After functional minimization one obtains for the order 
parameter m = \JujYl < <f) > the differential equation 



Am = Tm + 2m 3 , 



(17) 



where m = m(r,9,T;^f). Since the scaling functions in 
Eq. |JI} have the limiting behavior shown in Eq. J5J, 
where (3 = v = 1/2 within the mean-field approximation, 
in order to derive a boundary condition for the numerical 
calculation of the order parameter close to the surfaces 
of the wedge, i.e., 9^0 for r fixed, we seek a solution 
for Eq. fTTjl in the form 

m(r,9,T) = A ^ r ' T > + B(r, r) + C(r, t)9 + D(r, r)9 2 



(18) 



for 9 <C 1 (suppressing the 7-dependence in the notation) 
and obtain for both r > and r < 



1 n ^ 11 ( 1 Tr 



(B(r,r) = 0, D(r,r) = 0). This result agrees with the 
direct expansion of the mean-field profile obtained for 
r = [20]. In terms of the scaling functions P± this 
implies [see Eq. |JSJ|]: 



P + (r + ,0->0) = c+ 



P_(r_,6» -> 0) = c_ 



1 



r+ sin 9 



1 - -rlt 
6 + 



1 r 2 

12 " 



e>(6> 4 ) 



(20a) 



0{6 4 



+ O{0 3 ) (19) 



(20b) 

with c + = v/2 and c_ = 2 @. 

We use a numerical method |35l l3^ | to minimize Eq. 
(I16|l with respect to the order parameter profile at a fixed 
temperature, which is then subsequently varied. For 
computational purposes we choose suitably shaped finite 
volumes V w ( r ) for different opening angles 7 of the wedge 
(w) or ridge (r). We refrain from describing this choice 
of volumes here, because it does not matter in the calcu- 
lation of the universal amplitude functions T±(7). The 
choice of the volume is relevant only for that line con- 
tribution to the excess adsorption that depends only on 
the order parameter profile close to an infinite planar wall 
and thus is independent of the opening angle of the wedge 
(see the Appendix). As the temperature is changed we 
rescale the size of the volume V w ^ in accordance with 
the change of the correlation length £ = £ot~". This way 
we control the finite size effects caused by the finiteness of 
V w f r ) dictated by computational necessity. The finite size 
effects manifest themselves even close to those boundaries 
of the chosen volume that are far away from the walls of 
the wedge because of using approximate boundary con- 
ditions at these boundaries (see below). By effectively 
increasing the rescaled volume upon approaching T c , the 
values of the profiles at fixed spatial points within this 
rescaled volume converge to a limiting value. 

We choose a two-dimensional grid (i, j) and calculate 
the deviation of the order parameter profile from the 
known profile at T = T c at the given grid points. The 
profile near the confining surfaces of the wedge (ridge) is 
fixed according to Eq. (|19|) . (The grid points are lined 
up parallel to the walls of the wedge.) At the surfaces 
of the finite volume V w r r \ that are further away from the 
walls of the wedge, we prescribe initial values of the pro- 
file. Keeping these values fixed we then calculate new 
values of the profile inside the volume using the method 
of steepest descent. Having obtained the new values for 
the profile close to these surfaces, we change the profile 
at these surfaces proportionally to the change in their 
neighborhood, if it changes significantly in the direction 
perpendicular to the surfaces (the general case), or set it 
equal to that in neighboring layers if the profile is approx- 
imately constant close to these surfaces (for example far 
away from the edges of the wedges parallel to the con- 
fining walls of the wedge). The rules for the iteration 
according to the method of steepest descent in the space 
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of the parameters 

ay = m{rij , % , 0) - m(rij , 0y , t) (21) 



are 



fn+l) 



» 



dH(c 



da; 



(22) 



where k is a convergence parameter. With this method 
it is not necessary to calculate the Hamiltonian itself but 
only its derivative. The method has been described in 
detail in Refs. [3f| and (3^, so here we only want to 
point out that the divergence of the profile close to the 
walls of the wedge causes a divergence in the derivative 
of the Hamiltonian, too. This can be avoided, if one 
takes into account the known form of the divergence of 
the profile close to the edges of the wedge [see Eq. <|19[) ]. 
To this end we write ay as a product of two terms, one 
of which we choose to be such that when multiplied by 
the profile close to the boundaries (in the calculation of 
the derivative) it cancels the divergences of the profiles 
[see Eq. (|19|) ]. yielding a smooth gradient in the parame- 
ter space. We approximate the integrand of the gradient 
of H{cj}} [see Ea. lfTB|) ] by a sum of delta- functions posi- 
tioned at the grid points, so that the integral reduces to 
a simple sum over these points; for each opening angle k, 
is optimized separately for best convergence. 

The scaling function of the order parameter profile in 
a wedge [see Eqs. Oh©] with an opening angle of 90° is 
shown in Figs. 12151 One can easily see that the contour 




FIG. 2: The scaling function P+ of the order parameter profile 
in a wedge of opening angle 7 = 7r/2 with the edge located at 
the back of the figure at x' = f = {x' = rcos#, £ = rsinO 
[see Figs H1and ll4f all 1. The positions of the walls of the wedge 
coincide with the coordinate axes as indicated here with broad 
lines. The values of P+ at which contour lines are drawn are 
multiples of 0.1 and range from 0.1 to 1.0. 

lines quickly become parallel to the walls of the wedge 
as we move away from the edge (see Figs. |21 and 121 . 
This is especially apparent as we approach the walls. As 
one moves along the bisector of the wedge, the maximal 
curvature k of the contour lines decreases sharply (Fig. 
|3J). This underscores that in terms of the rescaled vari- 
ables, to a good approximation the effects of the edges 



are spatially localized. The maximal curvature k depends 
linearly on the values of P+ corresponding to the contour 
lines within the range of P+ values analyzed in Fig. [3] 



+ 




X 



+ 



FIG. 3: Projection of the contour lines of the scaling function 
P+ of the order parameter profile in a wedge of opening angle 
7 = 7r/2 onto the (x' + X+) plane. The edge is located at 
x' + — = 0. The curves correspond to values of P+ ranging 
from 0.1 to 0.6 with an increment of 0.1 from top to bottom. 
The inset shows that the maximal curvature n of the contour 
lines, occurring on the bisector (•), depends linearly on the 
values of P+ corresponding to the contour lines within the 
range of P+ values considered here. 

Along radial directions, i.e., for = const. (Fig. 0J) 
the scaling function exhibits a power law limiting be- 
havior close to the walls in accordance with Eq. (JSJ) 
and the paragraph preceding it. For large r± the be- 
havior crosses over into an exponential decay: P + (r + — > 
oo,0; 7) = A + (9 7 j)e- r + sint> , P_(r_ -> 00, 0; 7) - 1 = 
A_(0, r y)e~ r ~ Slne , where near the walls A± reduce to the 
amplitudes of the exponential decay away from an infinite 
planar wall: A+(6 -> 0,7) = V8 and A-{6 ->■ 0,7) = 2. 
The dependence of A±(9, 7) on is weak. The latter val- 
ues are valid for < 30°, beyond which the prefactors 
of the exponential functions increase with 0. Upon ap- 
proaching the walls of the wedge vertically, i.e., for 6 — > 
or 7 with r± fixed [see Eq. JHJ) and Fig. |S] the divergence 
of the profile has a power law form. 

The limiting behaviors of the scaling functions P± close 
to the edge of the wedge (r± — > 0) are described by the 
amplitude functions c±(0, 7) [see the paragraph following 
Eq. Q and Fig. EJ: 



P±(r±^ 0,0) = c±(0, 7)7-7^. 



(23) 



These functions are plotted for 7 = 7r/2 in Fig. H3 Ac- 
cording to Eq. JHJ, close to the walls of the wedge, i.e., 
for 0^0 these functions are given by c±(0 — > 0, 7) = 
c ± (sin0)-' 3 / l/ . 

The scaling function of the order parameter profile at 
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FIG. 4: The scaling function P+ of the order parameter as 
a function of the distance r+ from the edge of a wedge of 
opening angle 7 = tt/2. The curves correspond to values of 9 
ranging from 5° to 25° with an increment of 5° from top to 
bottom. 
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FIG. 5: The scaling function P+ of the order parameter along 
arcs of different radii r+ centered at the edge of a wedge of 
opening angle 7 = 7r/2. The curves correspond to values of 
r+ from 1 to 5 with increments of 1 from top to bottom. The 
curves are symmetric around 6 — 7/2 = 7r/4 and diverge as 
(c+r^^ l ')6~ l3 / v for 6 — » 0. For r+ = 1 the comparison with 
the asymptotic behavior c+(r+ sin Q)~^^ v for 6 < 7/2 and 
c+[r+ sin(7 — 9)]~ 13 ^ 1 ' for 6 > 7/2 is shown as a dashed curve 
[see Eq. ©]. 



a ridge [see Eqs. Q)-©] with an opening angle of 260° 
is shown in Figs. [7| and |SJ One can easily see that the 
contour lines rapidly become parallel to the walls of the 
ridge as one moves along them further away from the 
edge. As one moves along the bisector of the ridge away 
from the edge, the maximal curvature k of the contour 
lines decreases (Fig. [SJl. n depends linearly on the values 
P± characterizing the contour lines and thus has similar 
limiting behaviors close to the edge of the ridge and far 
from it as a function of r±, at least within the range 
studied in the inset of Fig. |SJ 

Comparing the order parameter profiles in a wedge or 




FIG. 6: The amplitude functions £±(#,7 



[P±(r±,t 



TT/2) 



in a wedge of opening angle 7 (t > 



: +, t < : x). The solid line corresponds to the func- 
tion c+ismey 13 ^ (t > 0) [see Eq. and the dotted line 
to c_ (sin ff)~^l v (t < 0), which are valid in the asymptotic 
regime 9 — > 0, but provide a surprisingly good description 
throughout the whole angle range < 6 < 7r/4 ~ 0.785. 
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FIG. 7: The scaling function P+ of the order parameter profile 
at a ridge of opening angle 7 = 260° with the edge perpen- 
dicular to the (x+, z+) plane and located at x+ — z+ = [see 
Figs. Q and, c.f.. I14f b11. The positions of the walls of the 
ridge are indicated here with broad lines. The values of P+ at 
which contour lines are drawn are multiples of 0.1 and range 
from 1.0 to 0.1 (top to bottom). 



at a ridge with the profile near a planar wall, one can 
visualize the wedge or ridge as being formed by breaking 
the planar wall into two halves, which in the case of a 
wedge are brought closer to each other, and in the case 
of the ridge are taken further apart. Close to the edge of 
the wedge this increases the values of the profiles, while 
close to the edge of a ridge these values are decreased as 
compared to the profile near a planar wall. 



B. Excess adsorption 

The presentation of the full order parameter distribu- 
tion requires to keep track of four variables: r,8,t, and 7. 
Therefore it is advantageous to consider also the excess 
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x + 

FIG. 8: Projection of the contour lines of the scaling function 
P+ of the order parameter profile in a ridge of opening angle 
7 = 260° onto the (x+,z+) plane. The projection is shown 
only on one side of the bisector of the ridge, which coincides 
here with the z+ axis. The wall of the ridge is indicated with 
a broad line. The edge is located at x+ — z+ = 0. The curves 
correspond to P+ = 0.1, 0.5, and 1.0 from top to bottom. The 
inset shows that the maximal curvature k of the contour lines 
occurring on the bisector is a linear function of P+ within the 
range considered here. 



adsorption in wedges and at ridges, which is experimen- 
tally relevant and provides reduced information depend- 
ing only on t and 7 [Eq. l|l(J|l ]. We are particularly in- 
terested in the line contribution [see Eqs. i|ll|) and Q13JI] 
characterizing the effect of the wedge (ridge) geometry 
via its universal amplitude functions [Eq. (|15J) ] . In order 
to calculate this quantity we use the fact that the integral 
in Eq. I|15|l can be rewritten as 

/>OG 

r±(7) = / F(P ± )dP±, (24) 
Jo 

where F(P±) is the area enclosed by the contour lines of 

P± = P±(r±,6;-y)-P± /2 (C±(r±,6)) (see Fig. £). Based 
on these areas one is left with a one-dimensional integra- 
tion to obtain T± numerically. Furthermore, exploiting 
the observation that the geometrical shapes formed by 
the contour lines are similar to one another for small and 
large areas, respectively, and using the limiting behavior 
of the scaling functions, we approximate F(P±) for small 
values of P± in terms of powers of P± , and for large P± 
in term of powers of ln(P-t ) . These approximate power 
laws are calculated based on different intervals in P± 
chosen as ever narrowing slices of that interval in P± in 
which the numerical data lie. The narrowing intervals ap- 
proach the small and large P± limit, respectively. With 
these power law approximations for different intervals in 
P± we obtain a series of approximate integrals for those 
intervals, for which due to practical limitations there are 




FIG. 9: Contour levels of the integrand of the integral that 
determines the amplitude of the reduced excess adsorption for 
T > T c [see Eqs. I|15fl an d EH hr a wedge of opening angle 
7 = tt/2 with the edge located at the back of the figure at 
x' = ( — (x' — r cos 6, £ = r sin 9). The walls of the wedge, 
indicated by broad lines, coincide with the coordinate axes. 
The values of P+ at which contour lines are drawn increase 
with multiples of 0.1 and range from 0.05 to 1.05. 

no numerical data (for small and large values of P±), and 
take the limit. This enables us to carry out the integral 
in Eq. H'24|l numerically for the whole range from P± = 
to P± = 00. 

This integration leads to the universal amplitudes 
r±(7) as shown in Fig. ^3 For small opening angles 




0.5 1 1.5 2 

y/tt 

FIG. 10: The amplitude of the line contribution to the excess 
adsorption for t > [Eqs. ill 1 fi and Ijl5^ ] as a function of the 
opening angle 7. The full line is a fit through the data (+). 
The inset shows that over a wide range F+(7) ~ 1.62/(7/77). 
F + (tt) = by definition. Note that T + (y) < |r+(27r - 7)|. 

7 they diverge as T±( r y — > 0) ~ I/7, vary linearly close 
to 7 = 7T for 7 < 7T, and their absolute values increase 
rapidly for 7 > tt. Numerical evidence suggests that this 
latter increase is exponential (but no divergence [see Fig. 

cm. 

Strikingly, the reduced excess adsorption r+(7) above 
the critical temperature appears to be proportional to 
the reduced excess adsorption below the critical 

temperature. We have calculated their ratio for seven 



9 



opening angles ranging from 20° to 240° and found 

r+(7)/r_(7) = 1.137 ± 0.006, (25) 
i.e., their ratio appears to be independent of 7. 



Ratios of the amplitudes of the excess adsorption above 
and below T c have been investigated for the case of a 
planar wall theoretically 0, experimentally Q, and us- 
ing Monte Carlo simulations The values obtained 
experimentally for the ratios of the amplitudes for the 
planar case (mean value: 1.19 ± 0.04) agree rather well 
with the result of the Monte Carlo simulations (1.11); the 
corresponding mean field value is l/y/2 @- 
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The angular dependence of the reduced excess adsorp- 
tion shows that for large angles, i.e., for a ridge, the ab- 
solute values are larger than for small angles, i.e., for 
a wedge: r±(7) < |r±(27r — 7)|. One can calculate 
the excess adsorption at a periodic array of wedges and 
ridges (see, c.f., Fig. 113(1 if their edges are sufficiently 
far apart from each other so that their influences do 
not interfere (see part B of the Appendix). The result 
expressed by Eq. IA.5I in the Appendix shows that in 
this limiting case the line contribution to the excess ad- 
sorption, which captures the effect of the wedges and 
ridges, is the sum of the contributions of single wedges 
of opening angle ^ w and of single ridges of opening an- 
gle 7 r = 2ir — j w sharing the same temperature depen- 
dence |t| /3_2ly . Thus the amplitude r ± of the combined 
contribution of one wedge and a neighboring ridge form- 
ing the basic building block of the array is given as as 
r ± = T±("/ w ) + T±("/ r ). This quantity may be viewed 
as a function of j r — -f w = 2ir — 2j w = 2(n — 7^), which 
characterizes the roughness of the surface. r + (7,. — j w ) 
is plotted in Fig. ^2 for (t > 0). One can see that all 
values are negative, i.e., the total excess adsorption (rel- 
ative to a planar substrate with the same area as the 
actual one of the corrugated surface) is decreased by the 
line contribution. This demonstrates that the decrease 
in adsorption for a ridge with j r = 2tt — j w dominates 
the increase due to the corresponding wedge with open- 
ing angle 7^ (see Fig. I10|l . The amplitude r + varies 
quadratically for small roughness and exponentially for 
large roughness. 



In Ref. [23 the total excess adsorption has been calcu- 
lated for curved surfaces. For a curved membrane with 
both sides exposed to a fluid near criticality, the sum 
of the excess adsorptions on the two sides per unit area 
was found to be larger for spherical regions of the mem- 
brane and smaller for cylindrical regions as compared to 
that for flat regions. Since the cylindrical regions may be 
viewed as rounded wedges and ridges, this latter finding 
exhibits the same qualitative trend as found here for the 
periodic array of wedges and ridges. 



FIG. 11: The amplitude F + of the line contribution to the 
excess adsorption of the basic unit of one wedge of opening 
angle y w and one ridge of opening angle 7 r = 2n — 7™ in 
an array of effectively independent wedges and ridges as a 
function of (y r — y w )/n, which characterizes the roughness 
of the surface (see, c.f., Fig. The dashed lines are fits. 

The insets show the behaviors in the limiting cases on a log- 
log scale and a log-linear scale for small and large roughness, 
respectively. They indicate that T + ((7,. — j w ) — > 0) ~ (7,. — 
7„) 2 and that F™ r increases as exp{[(-y r — lw)/^] 3 } up to 

(7r - 7w)/7T = 2. 



IV. FREE ENERGY OF THE CONFINED 
FLUID AND CASIMIR TORQUE 



In the previous chapters we have investigated struc- 
tural properties, i.e., order parameter profiles of critical 
fluids in wedges and close to ridges. In the followings we 
comment on their thermodynamic properties based on 
the free energy of such systems. 

The volume V(j) of the system shown in, c.f., Fig. ^| 
(for another possible choice see, c.f., Fig. I14f> is bounded 
by the walls of the wedge or ridge, which in turn end in 
the linear edge forming the vertex. Accordingly in the 
thermodynamic limit the total free energy F decomposes 
into a bulk, a surface, and a line contribution: 

F(T; 7 ) = V( 7 )/ 6 (T) + Sf s (T) + Lf t (T, 7), (26) 

where fb(T) is the bulk free energy density, f s {T) is the 
surface free energy density, and fi(T, 7) is the line free 
energy density. S is the total surface area of the wall in 
contact with the fluid, L is the length of the edge. Each of 
the three terms in the free energy and thus the total free 

energy itself are sums of a singular part (fl%ng\t,'y)), 
which contains the thermodynamic singularities in the 
vicinity of the bulk critical point t = (T — T c )/T c — > 0, 

and an analytic background contribution (/^'l' (T, 7)). 
The leading singular part of the total free energy can 
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be written in the form (see, e.g., Refs. 01 > an d E3) : 

|2-Q 





V(7) 




(C ± ) d 


S 












+ (tir 





(£o) d V a(l-a)(2-a) 



\t\- 



a s (l - a«)(2 - a s ) 



of (7) 



a;(l - a/) (2 - a z ) 



1*1 



2-a. 



(27) 



Here a ~ 0.11 is the bulk specific heat exponent, a s = 
a + v, ai = a + 2f , and af and are universal bulk 
and surface amplitudes. The bulk contribution depends 
trivially on 7 via the geometry ^(7), whereas the surface 
contribution is independent of 7. The line contribution 
carries a nontrivial dependence on 7 via the universal 
amplitude functions 0^(7). 



The background contribution takes on the form 



"back(T, 7) 

k B T c 



s 



1 fback (T) 



(28) 



If one of the sidewalls is moveable around the vertex with 
the far end suspended at, say, a force microscope, the 
torque 



M = - 



dF(T; 7 ) 
d-f 



(29) 



exerted by the fluid in the wedge or ridge on its sidewalls 
is experimentally accessible: 



M ± 



\tr a +fSL(T) 



a(l - a) (2 - a) 



-(d-2) 



daf(j)/d'y 



\2-a, 



a{){2 
I 



m) 



(30) 



With the bulk contribution known independently, this 
measurement provides access to the universal amplitude 
functions 0^(7) by focusing on the thermal singularity 
~ |t| 2 ~ Qi = \t\ v = \t\ - 63 in d = 3, since dv = 2 - a. 
The singular contribution to M can be called a critical 
Casimir torque. For fluids information about the back- 
ground term can be obtained from Eq. (A7) in Ref. |4S[ 
and from Refs. (HHSIp. 

Within mean-field theory 2 — ui = 1 so that the singu- 
lar line contribution becomes indistinguishable from the 
analytical background contribution. Thus our present 
approach renders only the sum of these two contributions 
without the possibility to isolate the amplitude functions 
Oj (7). As indicated by the pole ~ 1 } a for d— > 4, inclu- 
sion of Gaussian fluctuations beyond the simple mean- 
field theory is expected to generate a term ~ t In \t\ due 
to the resonance of the singular contribution ~ ■*-= — t 2 ~ a ' 
with an analytical background term ~ t [521]. The am- 
plitude of the singular term ~ iln|i| would allow one 
to retrieve at least the mean- field expression for 0^(7). 
However, this technically challenging inclusion of Gaus- 
sian fluctuations is beyond the scope of the present work. 

A suitable approach to obtain the change of the free 
energy upon varying the opening angle of the wedge or 
ridge involves calculating the field theoretical stress ten- 
sor. Analogously to the free energy density, the corre- 
sponding torque requires additive renormalization up to 
second order in temperature ■ We have followed this 
route within the present mean-field theory without iso- 



lating the critical Casimir torque. We have found that 
this combined torque diverges as M ~ I/7 2 for small 7, 
and it appears to be a linear function of I/7 2 within the 
angle range between 7 = 20° and 7 = 280°. 



V. SUMMARY 

In the present study of critical adsorption in wedges 
and close to ridges (see Fig. we have obtained the 
following main results: 

(1) We have discussed the scaling properties of the or- 
der parameter profile m±(r 1 1; 7) = a\t\P P±(r /£±, 6; 7) in 
terms of the bulk correlation length £± = above 
and below the critical point T c with t — (T — T c )/T c . 
The universal scaling functions P±(r/£±, 9; 7) diverge ac- 
cording to a power law close to the walls of the wedge or 
ridge, and decay exponentially far away from the walls 
[Eq. ©]. 

(2) In the thermodynamic limit the excess adsorption 
r±(sxj sh , t; 7) [Eq. (|10|) ] for volumes with linear exten- 
sions s± and S|| (see Figs. 1121 and m )l decomposes into 
a surface contribution that scales with the actual surface 
area of the confining walls (s_i_sjj 2 ') and a line contri- 
bution that scales with the extension in the invariant 
directions (sjj 2 ^) [see Eq. JUJ] as described in detail 
in the Appendix. The line contribution is the specific 
contribution arising from the influence of the edge on the 
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order parameter profile. Its amplitude r ; (t, 7) has the 
scaling form Tf(t,j) = a£^ 2 |t|^- 2l T ± (7) with the uni- 
versal amplitude functions T±(j) [Eq. 1)15(1] carrying the 
dependence on the opening angle 7. 

(3) We calculate the above scaling functions within 
mean-field theory [see Eqs. 1(16(1 and JT7J|] using a numer- 
ical algorithm both above and below the critical temper- 
ature (for a wedge see Figs. 12151 for a ridge see Figs. 
and|BJ. The amplitude functions of the power law diver- 
gence of the profile close to the walls [Eq. ©] are shown 
in Fig. El 

(4) Our numerical calculation also yields the experi- 
mentally relevant excess adsorption within mean-field ap- 
proximation. The universal amplitudes r + f/y) are shown 
in Fig. I1UI For small opening angles 7 they diverge as 
r + (7 —> 0) ~ I/7, vary linearly close to 7 = n for 7 < it, 
and their absolute values increase rapidly for 7 > tt. Nu- 
merical evidence suggests that this latter increase is ex- 
ponential, but without divergence [see Fig. ITT) . The 
reduced excess adsorption T + (7) above the critical tem- 
perature appears to be proportional to the reduced excess 
adsorption T_ (7) below the critical temperature with 
r+(7)/r_(7) = 1.137 ± 0.006. We have considered a 
wedge and a ridge together as forming the basic unit in a 
periodic array (see Fig. I13fl . The total excess adsorption 
relative to that of a planar substrate with the same area 
as the actual one of the corrugated surface [Eq. (|A.6ll ] is 
decreased by the line contribution (see Fig. 111(1 . 

(5) The variation of the free energy of the system with 
the opening angle of the wedge or ridge gives rise to a 
torque acting on the sidewalls [Eq. (J21J)] . The free energy 
decomposes into a singular contribution exhibiting scal- 
ing [sec Eq. (|27|) ]. and an analytic background contribu- 
tion [Eq. i(28|) ]. In d = 3 the critical Casimir torque varies 
as a ; (7)|t|" with universal amplitude functions 0^(7). 
This cusplike temperature singularity is expected to be 
experimentally accessible via suitable force microscopy. 
The theoretical calculation of the corresponding ampli- 
tude functions remains as a challenge. 
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APPENDIX: DECOMPOSITION OF THE 
EXCESS ADSORPTION 

In this appendix we discuss how the excess adsorption 
in a wedge or at a ridge decomposes into surface and line 
contributions [see Eqs. (|TT|) - ((T3*|) ]. In the definition of 
the excess adsorption [Eq. 110(1 ] one considers a finite 
volume V of integration that is enlarged to fill the total 
volume of the wedge or ridge in the thermodynamic limit. 
As shown below for two examples, the expression for the 



line term (such as Eq. ((13(1 ) actually depends on the 
choice of the shape of the volume V. 

In the following we first analyze the excess adsorption 
for a single macroscopic wedge or ridge which will 
be followed by a discussion of the excess adsorption for 
two possible experimental realizations (|"2l and H?)) . 



1. A single macroscopic wedge or ridge 

a. First choice of the volume of integration 

Our first choice for V is shown in Fig. ^|for the case of 
a wedge. This choice is inspired by the idea that the sin- 
gle wedge or ridge considered here is ultimately a member 
of a periodic array. With this polygonal cross section and 
a similarly constructed one for the ridge, one can cover 
the total volume Vtot of a fluid in contact with a sur- 
face formed as a periodic array of wedges and ridges in 
a natural way (see Fig. 113(1 . All the formulas explicitly 
stated below for the wedge are valid for the ridge, too. 
The volume V is symmetric with respect to the biscc- 



Z Jl 




FIG. 12: One choice for the volume V of integration in the 
definition of the excess adsorption [Eq. HUI l with its cross 
section in the plane perpendicular to the invariant directions. 




FIG. 13: A fluid is bounded by a substrate shaped as a peri- 
odic array of wedges and ridges with a finite depth D. The 
total volume of the fluid Vtot can be naturally decomposed 
into subvolumes of the type shown in Fig. 1121 

tor plane of the wedge, which allows us in the following 
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to consider only one half of the wedge, and multiply the 
corresponding expression by a factor of two. 

As indicated in Fig. El the volume V in Eq. itTTl 
is finite and thus the integral is finite, too. As a first 
step in carrying out the thermodynamic limit we keep 
the shape of V but shift the upper boundary z = const 
to infinity. Since m± approaches mi, exponentially this 
extension of V increases r± by an exponentially small 
amount and thus does not contribute to the two leading 
terms under consideration in Eq. (|llfl . (In the spirit 



of the thermodynamic limit one first increases V before 
one can possibly consider the limit t — > 0. Therefore 
these arguments are not impaired by a power law decay 
of m±(z — > oo,t = 0) in the thermodynamic limit.) For 
the resulting semi-infinite strip Eq. (|10f) can be rewritten 
in the following form by adding and subtracting the order 

parameter profile mj. (C = r sm ®i t) of a fluid in contact 
with an infinite planar wall in the integrand (see Figs. ^ 
and ITU): 



f±(fl X ,*i|,t;7) = 2s t 2 



d( 



2 + tan( 7 /2) 



c 



dx'[m ± (r(x', O,0(x', Q, t; 7) - m± /2 (C, *)] 



-2s 



d-2 



tan( 7 /2) 
- + 



<K 



— -t-tan( T /2) 



tan(-y/2) 



dx'[m¥ /2 ((,t)-m b (t)} 



(A.l) 



where x' = rcos8 = z cos(7/2) + x sin(7/2) is the coordi- 
nate measured from the apex parallel to the nearest wall 
of the wedge and ( = rs'm8 = zsin(7/2) — xcos(7/2) is 
the normal distance from the nearest wall of the wedge 
(see Fig. EJl. 

In the inner integral of the first term the upper inte- 
gration limit can be shifted to infinity, i.e., sj_ — ► +oo, 
with an addition of only exponentially small corrections 
to the integral, because m± approaches m± ^ 2 exponen- 
tially for x' — > +oo at fixed £. Thus the first term in Eq. 
(IA.1II approaches a constant for s± — ► +oo and this con- 
stant involves an unlimited integral over the whole half 
of the wedge. Expressed in terms of cylindrical coordi- 
nates this term yields the line contribution T^(t, j)s^ 2 ^ 
in Eqs. (fTTfl and iQ . 

As the integrand of the second term in Eq. (|A.1() does 
not depend on x', the inner integration simply yields a 
factor sj_/2. The outer integral then yields in Eq. 

©• Together with Eq. (JT2) this verifies Eq. (fTT|) . 



b. Second choice of the volume of integration 

Our second choice of the shape of the volume V as 
shown in Fig. I14f a1 corresponds to the one used in Ref. 
|37j . where liquids confined by two opposing structured 
walls have been studied; in this geometry one cannot 
infinitely extend the volume in the +z direction. On 
the other hand in the case of the ridge as shown in Fig. 
I14f b). starting out from a finite volume V, the upper 
boundary of the volume [z = const) can be shifted to 
+oo as well as the two vertical boundaries (x = ±const) 
to ±oo with only exponentially small corrections to the 
integral in Eq. (JTUJ, because in the directions z — > oo 
and x — > ±oo for z fixed the order parameter attains 




FIG. 14: A second choice for the volume V of integration 
in the definition of the excess adsorption [Eq. 1KH ] for a 
wedge (a) and a ridge (b) with its cross section in the plane 
perpendicular to the invariant directions. 



its bulk value exponentially. As for the previous choice 
also for the present geometries Eq. 1|10|) can be rewritten 
by adding and subtracting the order parameter profile 
m«^ 2 (C _ rs [ n Q f^ f a fl^d m contact with an infinite 
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planar wall in the integrand (see Fig. I14fl : 



T±(s ± ,s h t;~f) = 2s 



d-2 



-t sln 7 



<K 



+2s 



d-2 



~t sin 7 



^-Ctan( 7 /2) 



dx'[m±(r(x',(), 8(x', Q,t; 7) - m± /2 (C, *)] 



-Ctan( 7 /2) 



rfx'[m^ /2 (C,i)-m b (f)] 



(A.2) 



where as before x' = r cos is the coordinate parallel to 
the wall, and C, = rsin# normal to it [see Fig. I14f a)]. 
In the first term the upper integration limits of both in- 
tegrals can be shifted to infinity, i.e., s± — ► oo, with an 
addition of only exponentially small corrections to the 
integral, because m± approaches exponentially for 

x' — > +oo at fixed and m± attains its bulk value ex- 
ponentially for £ — > +oo at fixed x' . Thus for the first 
term the limit sj_ — > oo exists and is finite with the two- 
dimensional integral covering the whole half of the wedge. 
Expressed in terms of cylindrical coordinates this term 



with 



yields the line term T^(t,j)s 



[see Eqs. JTIJ and 
(II 31) ] as for the previous choice of the volume. 

The inner integral in the second term can be car- 
ried out, because the integrand is independent of x'. 



c 



tan 5j- 



tan 7/2 



This yields a prefactor ^ 

i |sj_ — £j^C|- The first term of this prefactor gives 

rise to the surface contribution in Eqs. J2J, pijl. and 
(|12fl . if one shifts the upper integration limit of the £ in- 
tegration to infinity with the addition of exponentially 
small corrections. After multiplying this prefactor by 
two, its second term gives rise, however, to another line 
contribution with the amplitude 



sin 7 



([m^ /2 ((,t)~m b (t)}d<;, (A.3) 



where the upper limit of integration has also been shifted 
to infinity with an exponentially small correction. We 
note that this additional line term depends on the order 
parameter profile at a planar substrate only. Due to the 
extra factor £ in the integrand, the integral in Eq. l!A.3() 
is finite for d = 4 in spite of m^ 2 (C — * 0) ~ £ _1 m this 
case, i.e., T± does not carry a factor proportional to -^-^ 
as r[^ 2 does (compare Eq. JSJ). 

Thus in the thermodynamic limit the two choices of 
the integration volume V yield the same surface contri- 
butions to the excess adsorption but different sub- 
dominant line contributions Tf 1 : 



d d x m±(x) = Vm b + STf + LTf + ©(s^ 1 ) (A.4a) 



rf = r±(t, 7 ) = 



and 



p/2 roc 

2 d0 drr[m±(r,6,t;j) 
Jo Jo 

-m± /2 (({r,9),t)}, choice I, (A.4b) 



r± (*,7) = r±(t,7) + rf(*, 7 ), 



choice II, 

(A.4c) 

with S = s ± s|[~ 2 , L = s|[~ 2 , and ff given by Eq. (fO)> . 
Experiments cannot be carried out for infinitely deep 
wedges. Instead they can be carried out for either a 
periodic array of wedges or for a single wedge of finite 
depth carved out from a wide planar surface. In both 
cases additional ridges must be formed giving rise to their 
own adsorption properties. Therefore experiments on 
such systems give access only to certain combinations of 
wedge and ridge excess adsorptions, whose corresponding 
line contributions carry the relevant additional informa- 
tion about the nonplanar substrate geometry. The actual 
choice of the corresponding integration volume depends 
on the actual experimental setup (compare, e.g., Fig. 1131 
see Subsecs. I"2land r3*|) The results of this subsection show 
that the line contributions depend on such details even 
in the thermodynamic limit. 



2. Periodic array of wedges and ridges 

In this subsection we consider a substrate with a pe- 
riodic series of edges and wedges as depicted in Fig. ^] 
There is a variety of experimental techniques to create 
such kind of surface morphology. If the opening angle of 
the wedges is 7, the opening angle of the ridges is 2ir — 7. 
Here we focus on the limiting case that the depth of the 
wedges D = s±cos('y/2) is sufficiently large, so that the 
deviations of the profiles close to the edges of the wedges 
and ridges, respectively, from the profile of a fluid ex- 
posed to an infinite planar substrate do not influence 
each other. As apparent from Figs. 1121 and fHfl this case 
corresponds to the first choice of the volume of integra- 
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tion for the single wedge or ridge and thus leads to the 
following decomposition of the excess adsorption [see Eq. 
(jS3bJ] 

/ d d x m±(x, t; 7) = V tot m b (t) + S to tTf(t) 
Jv tot 

+L(N w Tf w {t, lw ) + N r T±(t, lr )) + Oisl 1 ) (A.5) 
I 



where 7^, is the opening angle of the wedge, and j r = 

2ir — T^,. This is in accordance with Eqs. (32)-(36) of 

Ref. |37j . where a different coordinate system was used 
0. 



3. A single wedge embedded into a planar wall 

In this subsection we consider a single wedge of opening 
angle 7 carved out of a planar surface, thus producing 
also two ridges of opening angles (3ir — r y)/2 (see Fig. HT)l . 
Due to the symmetry of the configuration we consider 




FIG. 15: A fluid is bounded by a substrate shaped as a wedge 
carved out of an infinite planar wall. The total volume of 
the fluid can be decomposed into the numbered subvolumes 
that emerge naturally due to subvolumes 3 and 5 used for 
describing a single ridge (first choice of the volume of integra- 
tion 1 1 al ) and a single wedge (second choice of the volume 
of integration (1 1 bl l. respectively. 

only half of the wedge and focus on the case that s± 
is sufficiently large so that the ridges and the wedge do 
not influence each other. One half of the total volume 
of the fluid is decomposed into six numbered subvolumes 
as shown in Fig. El We calculate the excess adsorption 
[Eq. Q1U|)] separately for each one of the subvolumes by 
adding and subtracting the profile of a fluid exposed to 
an infinite planar wall m±^ 2 in the integral of the order 
parameter as in Subsecs. ITlandrSI 



with the total surface of the substrate Stot — Ns±Su. 2 , 
where N = N w + N r is the number of segments of length 
s±j r±(t) = r± /2 (t) [see Eq. ©], L = sjj~ 2 , N w = N r 
are the numbers of wedges and ridges, respectively, and 



(A.6) 

I 

For the first volume there is no line contribution re- 
sulting from the integral over the difference of the actual 
profile m± from that in front of an infinite planar wall 

m±^ 2 , because this difference is exponentially small in 
the thermodynamic limit. There is also no surface term, 
because this first volume touches the substrate only at 
one point. However, following similar considerations as 
for a single wedge or ridge, this volume gives rise to a 
line contribution [see Eq. (|A.4a|) ] to the excess adsorp- 
tion [Eq. HI due to the deviation of m±^ 2 from nib 
with amplitudes: 

/>oo 

f± (t )7 ) = tan^p / C[™± /2 (C,*) - ro 6 (i)]d£ 

(A.7) 

where we have shifted the z — const boundary to +00 
with an exponentially small correction to the integral. 

In the thermodynamic limit (s±, s'j_ — > 00) the adsorp- 
tion profile m± in the second volume will tend exponen- 
tially to that of a fluid exposed to an infinite planar wall 
m^ 2 . Thus it supplies a surface term with the ampli- 
tudes Tf(t) = T± /2 (£) [see Eq. ©] for the area s^sjj" 2 
We note that this subvolume two when extended to infin- 
ity in the +z direction overlaps with the other half of the 
wedge, but this results in only an exponentially small 
correction to the excess adsorption, because the profile 
m± approaches the bulk value mj, exponentially with the 
distance £ from the wall. The subvolume two does not 
generate a line contribution. 

The third volume gives the same contributions [Eq. 
HA.4a|l ] to the excess adsorption as a single ridge of open- 
ing angle (3ir — 7) /2 with the first choice of volume of in- 
tegration shown in Fig. 1121 i.e., a surface term with the 
amplitudes Tf(t) = T± /2 (i) for an area s_Lsjf 2 [see Eq. 

ijSJ] and line term amplitudes T^(t,j r = (3tt — j)/2) [see 
Eq. (|A.6(l ] . Note that the volume of integration can be 
extended to infinity in the direction parallel to the bisec- 
tor of the ridge even though in this case there is an over- 
lap with the other half of the wedge, because this results 
in only an exponentially small correction to the excess 
adsorption as the profile m± approaches the bulk value 



f7„( r )/2 POO 

r it(r)(*.T»(r)) = 2 d9 drr[ml {r \r,9,t; lw{r) )-m^ /2 (ar,9),t)}, 

Jo Jo 
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m b exponentially in this direction further away from the 
ridge. 

Similarly to subvolume one, subvolume four does not 
generate a line term resulting from the integral over the 
difference of the actual profile m± from m^ 2 . There 
is also no surface term, but similar considerations as for 
a single wedge or ridge show that this subvolume gives 
rise to a line contribution to the excess adsorption with 
amplitudes of the form: 



lt 4 (*,7) 



tan 



7T — 7 



7 

tan — 
2 



( / ([m^ /2 ((,t)-m b (t)]dC, (A.8) 
'o 



once we have shifted the x = boundary to +oo with an 
exponentially small correction to the integral. 

The fifth subvolume together with its counterpart at 
x > gives the same contributions [Eq. (|A.4cj to the 
excess adsorption as a wedge of opening angle 7 with 
the second choice of the volume of integration shown in 
Fig. ll4T ah i.e., a surface term with amplitudes Tf(t) = 
j gec ^ c axea, s±s^~ 2 , and the line 

term amplitudes rj^^i, 7) from Eq. l|A.4cll . 

Finally the subvolume six yields only an exponentially 
small contribution to the excess adsorption. 

Adding up all contributions to the excess adsorption 
in this geometry one obtains 



d d xm±{x,t; 1 ) = V tot m b (t) + S tot Tf(t) + L 



Vtot 



Tfjt, lw ) + 2Tf(t, lr = (3tt - lw )/2) 



-4 {-^ - — - \ tan ^) ^ C[m± /2 (C, *) - m b (t)}d( 
\cos ^ sin7„ 2 2 / J 



(A.9) 



As in the case of a periodic array of wedges and ridges 
the line contribution to the excess adsorption contains 
a combination of wedge and ridge terms. Since different 
combinations thereof enter into the excess adsorption of a 
periodic array and of a single embedded wedge, measure- 
ments of both of them provide independent information. 
However, these configurations do not allow one to access 



the ridge and wedge contributions individually. In the 
case of a single embedded wedge the line contribution 
to the excess adsorption contains in addition the first 
moment of the order parameter profile of a semi-infinite 
planar system, which can be determined independently 
from the knowledge of m°°/ 2 (£). 
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